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Abstract. Pull-back transformations between Heun and Gauss hypergeo- 
metric equations give useful expressions of Heun functions in terms of better 
understood hypcrgeometric functions. This article classifies, up to Mobius au- 
tomorphisms, the coverings P 1 — > P 1 that yield pull-back transformations from 
hypergeometric to Heun equations with at least one free parameter (excluding 
the cases when the involved hypergeometric equation has cyclic or dihedral 
monodromy). In all, 61 parametric hypergeomctric-to-Hcun transformations 
are found, of maximal degree 12. Among them, 28 pull-backs arc compositions 
of smaller degree transformations between hypergeometric and Heun functions. 
The 61 transformations are realized by 48 different Belyi coverings (though 2 
coverings should be counted twice as their moduli field is quadratic). The 
same Belyi coverings appear in several other contexts. For example, 38 of the 
coverings appear in Herfutner's list of elliptic surfaces over P 1 with four sin- 
gular fibers, as their j'-invariants. In passing, we demonstrate an elegant way 
to show that there are no coverings P 1 — > P 1 with some branching patterns. 



1. Context and overview 
The Gauss hypergeometric equation 

and the Heun equation 

(1.2) d2y(x) | ( C d a + b-c-d+l \ dY(x) | abx-g Y(x) = 
dx 2 \x x—l x~t J dx x(x—l)(x — t) 

are canonical second-order Fuchsian differential equations on the Riemann sphere P 1 , 
with 3 and 4 regular singularities, respectively. Transformations among these equa- 
tions give identities between their standard hypergeometric and Heun solutions. 
For example, there is a single covering P 1 — > P 1 of degree 2 (up to Mobius trans- 
formations). It induces the classical quadratic transformations of hypergeometric 
functions, such as 



(1-3) 2 Fi 



2A, 2B 
A + B+\ 



1 ( A + ll+i AX(1 - X) 



Moreover, the same covering induces the well-known Heun-to-Heun quadratic trans- 
formation 17, Thm. 4.1], and an identification of the general 2~Fi{A, B; C | 4a;(l— x)) 
function with a standard local solution of Heun's equation with the parameters 
(t,q,a,b,c,d) = (™, 2AB, 2A, 2B, C, C) . These transformations are parametric, 
since they have at least one free parameter such as A, B. 
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The aim of this paper is classification of all parametric pull-back transformations 
between hypergeometric and Heun functions. The considered pull-back transfor- 
mations are of the form 

(1.4) «■—> y(z)^Y(x) = 6(x)y(<p(x)), 

where ip(x) is a rational function and 6(x) is a radical function, i.e., a product of 



powers of rational functions. Geometrically, transformation (1.4 1 lifts or pulls back 
a Fuchsian equation on the curve P* to one on the curve P*, along the covering 
ip: P* — > P*. The gauge pref actor 9(x) is usually chosen such that the pulled-back 
equation has fewer singularities and canonical values of some local exponents. 

Pull-back transformations between Gauss hypergeometric equations were re- 
cently classified by Vidunas [29] . Next to the classical quadratic, cubic and Goursat 
[5] transformations, a few sets of unpredicted transformations were found, including 
parametric transformations from hypergeometric equations with cyclic or dihedral 
monodromy. Moreover, the hypergeometric-to-Heun transformations without the 
prefactor 9(x) have been classified by Maier [TO]. In both classifications, the heart 
of the problem is determining the covering maps ip(x) that can appear. They are 
typically Belyi maps, in the sense that (apart from dull exceptions of Proposition 



2.3 here) they have at most 3 critical values on the Riemann sphere P*. In fact, the 
critical values of those fix) are typically the singular points z = 0, z = X, z = oo 
of the hypergeometric equation, and the branching points include the singularities 
x = 0, x = 1, x = oo (and x = t) of the pulled-back hypergeometric (or Heun) 
equation. The approaches of [El [29] include: 

(?) determining the branching patterns that (p can have; 

(ii) determining which of those patterns can be realized by a rational function 

(p(x); 

(Hi) normalizing the points x — 0, x = 1, x = oo of <p(x), and deriving identities 
between hypergeometric and Heun functions by identifying corresponding 
local solutions of thereby related differential equations. 

This article follows this strategy and the techniques of [29 to generate a complete 
list of coverings tp that can appear in parametric Heun-to-hypergeometric reduc- 
tions. We find 61 different transformations (excluding infinite families of pull-backs 
from hypergeometric equations with cyclic or dihedral monodromy [32]), realized 
by 48 different Belyi coverings. An explicit formula for each covering is given in 
Table [ij The Belyi maps are not normalized for Step (Hi). The induced identities 
between hypergeometric and Heun functions are thoroughly examined in the paral- 
lel article [3D]. Here we not concerned with the technical issues of determining the 
prefactor 0(x), identifying local solutions, symmetries of the hypergeometric and 
Heun equations, nor even introducing Heun functions. 

By the Grothendieck correspondence [22] any Belyi map (p : P* — > P* corresponds 
bijectively to a dessin d'enfant on P*, up to Mobius isomorphisms of both Riemann 
spheres P*, P*. Generally, the dessins are defined combinatorially as certain bicol- 
ored graphs. For our purposes, the dessins d'enfant of a Belyi map <p(x) is the graph 
on P* obtained as the pre-image of the line segment [0, 1] on Pj, up to isotopy. The 
vertices above z — are colored black, and the vertices above z = 1 are colored 
white. The order of each vertex is equal to the branching order at the corresponding 
x-point. Figure [T] depicts the dessins for all 48 encountered Belyi coverings. Most 
of the white points have order 2, and then they are not depicted. Black points of 
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Figure 1 . Dessins d 'enfant of the Belyi coverings for parametric 
Hcun-to-hypergeometric reductions 



order 3 or 4 are not depicted either, unless they are connected to a white point 
of order 1. A thin edge connects a pair of displayed black and white vertices. A 
thick edge connects two black points (either displayed or clearly branching) with 
an implicit white point somewhere in the middle. Each cell (i.e., a two-dimensional 
connected component of the complement on Pj. , possibly the outer one) represents 
a point above z — 00. The branching order of each cell is determined by counting 
the number of black points met while tracing a loop along its boundary. 

It is instructive to follow the branching orders and incidences on the dessins 
while following our classification of possible coverings in Tables [l}|3j In princi- 
ple, the pull-back Belyi coverings can be classified by generating and counting the 
dessins satisfying the suitable branching patterns. However, it is difficult to ensure 
completeness of a large list of dessins. We first computed the Belyi coverings explic- 
itly, then easily generated the required dessins by combinatorial consideration. For 
each possible branching pattern, there is at most one Belyi covering except for the 
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coverings H21 and H^. Therefore completeness and identification of the dessins is 
quickly established. The coverings if 2 i, #44 are defined over Q(y/— 3) and Q(i), 
respectively. All other coverings are defined over Q and R, hence their dessins have 
a reflection symmetry. The dessins for H21, H44 should actually be counted twice, 
as the complex conjugation gives non-isotopic dessins. The proper count of dessins 
and Belyi coverings is therefore 50, not 48. 

Many of the encountered Belyi coverings occur in other contexts, particularly in 
the theory of elliptic surfaces and Picard-Fuchs equations. The coverings from Hi 
to i?38 occur in Herfurtner's list [10] of elliptic surfaces with four singular fibers, up 
to Mobius transformations. The order of these coverings follows [Till Table 3], and 
the numbering is used in |19j where the corresponding pull-backs to Heun equations 
(specializable to Picard-Fuchs equations for the elliptic surfaces) are observed. The 
coverings Hi to H 6 have the maximal degree 12, and produce the Beauville list 
[2] of the coverings generating semi-stable elliptic surfaces with four singular fibers. 
Their branching orders above z = are all 3, and above z = 1 they are all 2, as can 
be seen from the dessins. The branching pattern of Hi is written by us as follows: 

(1.5) [2] 6 = [3] 4 = 9+1+1+1. 

The four singular fibers of the corresponding elliptic surface have the Kodaira types 
Ig,Ii,Ii,Ii. This covering is also described as a Davenport- Stothers triple [S3]: it 
can be written as F 3 /G 2 , where F, G are polynomials of degree 4 and 6 (respec- 
tively), such that the polynomial F 3 — G 2 has the minimal possible degree 3. 

A pull-back transformation defined over K can be nicely illustrated by subdivi- 
sions of the Schwarz quadrangle for the pulled-back Heun equation into Schwarz 
triangles for the initial hypergeometric equations, following [TT| 112] . In the hyper- 
bolic geometry setting, these are Coxeter decompositions [7J or divisible tilings [5] 
of a hyperbolic quadrangle into mutually similar hyperbolic triangles. We describe 
these picturesque illustrations in §4.3| and Figure [2j 

This article is structured as follows. Section [2] establishes pivotal lemmas on the 
behavior of singularities and local exponents of Fuchsian equations under pull-back 
transformations. Section [3] presents the main results in Tables [T[j4] and explains 
them (and the notation) in a few steps. Of the three mentioned generation steps 



(i)-(iii), the first step is elaborated in §^3.1 3.2, while computations for Step (ii) 



are reviewed in |4.1| Step (Hi) is thoroughly considered in the parallel paper [30] . 
Furthermore, prelates our classification to Herfurtner's list [TU] and Felikson's list 



of Coxeter decompositions [7J, and £4.4 examines the composite transformations. 
Section [5] presents an elegant approach to prove non-existence (or uniqueness) of 
Belyi coverings with some branching patterns, and applies it not only to the ob- 
tained list of branching patterns, but also to the Miranda-Persson classification [18] 
of K3 semi-stable elliptic surfaces with six singular fibers. 



2. Pull-backs and local exponents 
The singular points and the local exponents of Gauss hypergeometric equation 



(1.1 ) are usefully encoded in the Riemann P-symbol scheme 
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The local exponent differences at the 3 singular points are therefore 
(2.2) 1-c, c-a-b, a-b. 



Similarly, the Riemann scheme of the Heun equation (1.2) is 
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(2.3) 
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- d — a — b 
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The parameters a, b, c, d determine the local exponents, while the parameter q is 
accessory. In particular, the 4 exponent differences are 

(2.4) 1-c, 1-d, c + d- a-b, a-b. 

The Heun equation contains many interesting special cases, including the Lame 
equation [5]. The Heun equation and its solutions appear in problems of diffu- 
sion, wave propagation, heat and mass transfer, magneto-hydrodynamics, particle 
physics, and the cosmology of the very early universe. 

By E(a,/3,j) we denote a Gauss hypergeometric equation of the form (1.1) 
with the exponent differences (2.2) equal to a, (3, 7 in some order. Similarly, by 
HE(a, (3,j,6) we denote a Heun equation of the form (1.2) with its exponent dif- 
ferences equal to a, (3, 7, 5 in some order. These notations do not assign local 
exponents to particular singular points, nor they specify the accessory parameter q. 

The degree of a pull-back transformation (1.4) between Fuchsian equations is 
the degree of the rational function f{x). The existence of a pull-back from some 
E(ai, Pi, 71) to some HE(a 2 , (3 2 , 72, 5%) of degree D will be indicated by 

(2.5) E( ai ,Pim) ^- HE(a 2 ,f3 2 , l2 ,5 2 ). 

Sometimes the pull-back covering or the transformation will be indicated more 
specifically by a subscript on the degree D. Similarly, 

E{ot x ,Pi,^i) 4^- E(a 2 ,f3 2 ,j 2 ), HE(ai,f3 1 ,ji,S 2 ) ^ HE(a 2 , /3 2 , 72, 5 2 ) 

will indicate pull-back transformations between hypergeometric or between Heun 
equations. For brevity, we refer to these three types of transformations as Gauss-to- 
Heun, Gauss-to-Gauss (or just hypergeometric) and Heun-to-Heun pull-back trans- 
formations. In particular, the 3 quadratic transformations mentioned at the begin- 
ning of this article actually are: 

(2.6) £(1/2, a, /3) A E(a, a, 2/3), 

(2.7) #£(1/2, 1/2, a, /3) HE(a, a, (3, /3), 

(2.8) E(a, f3, 7) ^~ HE (a, a, 2/3, 2 7 ). 

As in the notation (a±, 71) (a 2 ,/3 2 , 72) of [29], the arrows follow the direction 
of the covering cp; — > P^. To emphasize: these notations indicate the existence 
of some differential equations with the stated exponent differences that are related 
by a pull-back transformation, rather than the existence of a pull-back between any 
equations with the specified exponent differences. 

Our classification is obtained by considering the behavior of singularities and 
local exponents of Fuchsian equations under pull-backs. Any transformation of the 
form ( 1.4 1 pulls-back a Fuchsian equation to a Fuchsian equation, usually with more 
singular points. To pull-back a hypergeometric equation to a Fuchsian equation 
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with just 4 singular points, special restrictions apply to the covering ip(x) and the 
hypergeometric equation. 

The following definitions are taken from |29j . An irrelevant singular point of 
a Fuchsian equation is a non-logarithmic singular point where the local exponent 
difference is equal to 1. For comparison, an ordinary (i.e., non-singular) point is a 
non- logarithmic point with the local exponents and 1, and an apparent singularity 
is a non- logarithmic singular point with the local exponents and an integer k > 1. 
A relevant singular point is one that is not irrelevant. Any irrelevant singular point 



can be turned into an ordinary point by a pull-back (1.4) which is prefactor-only, 
i.e., one with ip(x) = x. Hence, what is of primary importance is how many relevant 
singular points the pulled-back equation has. This number is affected only by the 
choice of covering f(x), and not by the choice of prefactor 9{x). 

The following two lemmas describe the crucial behavior of singularities and local 
exponents under pull-backs. 

Lemma 2.1. Let ip: P* — > V\ be a finite covering. Let E\ denote a Fuchsian 
equation on W\, and let E 2 denote the pull-back on Pi. of E\ by transformation 



(1.4). For any S £ ¥]., let k := ord v (P) denote the branching order of ip at S. 

(a) The exponents of E% at S equal kct\ +7, ka 2 + 7, where: 

o a\,a 2 are the exponents of E\ at (p(S) G ¥],; 

o 7 is the exponent of the radical function 0(x) at S . 

(b) If tp(S) is an ordinary point of E\, then S will fail to be a relevant singular 
point for E2 if and only if k = 1 (i.e., the covering ip does not branch at 
S, i.e., S is not a branching point of ip). 

(c) If <p(S) is a singular point of E\, then S will fail to be a relevant singular 
point of E 2 if and only if 

o k > 1 and the exponent difference at <p(S) is equal to 1/k; or, 
o k = 1 and <p(S) is irrelevant. 
In either case S will be an irrelevant singular point or an ordinary point. 

Proof. The first statement is mentioned in the proof of j29j Lemma 2.4]. The other 
two statements are parts 2 and 3 of Lemma 2.4]. □ 

Lemma 2.2. Let ip: P^, — > V\ be a covering of degree D, and let A denote a set 
of 3 points on V\ . 

(a) If all branching points of (p lie above A, i.e., no point of P^\ A is a critical 
value of ip, then there are exactly D + 2 distinct points on P;j. above A. 
Otherwise, there are more than D + 2 distinct points above A. 

(&) If there are exactly D+3 distinct points above A, there is only one branching 
point that is not above A. 

Proof. The first statement is part 1 of [29, Lemma 2.5]. It follows from the Hurwitz 
formula [21 Corollary IV. 2. 4], which says that the sum of ord v (P) — 1 over the 
branching points f £ P], must equal 2(D — 1). The second statement is a slight 
extension (utilized in |13j). □ 

Suppose one starts with a hypergeometric equation E\ on P^ . Let A denote the 
set {0, 1, 00} containing the singularities of E\. It follows from the above lemmas 
that to minimize the number of singular points of a pull-back of Ex, one should 
typically allow branching points of ip only above A. Otherwise, there would be 
more than D + 2 distinct points above A, and generically, each of these D + 2 
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points would be a singular point of the pulled-back equation. By Lemma |2.1[ c), 
further minimization is possible if one or more of the exponent differences of E\ 
in A are restricted to be of the form 1/fc. 

Recall that a covering <p : P 1 — > P 1 is a Belyi covering [24] if it is unbranched 
above the complement of a set of three points, such as {0, l,oo}. By the above 
consideration, one expects that the pull-back coverings for Gauss-to-Heun transfor- 
mations will typically be Belyi coverings. The following proposition classifies the 
rather degenerate situations in which non-Bclyi coverings can occur. 



Proposition 2.3. Suppose there is a pull-back transformation (1.4 1 of a hypergeo- 
metric equation E\ to a Fuchsian equation with at most 4 singular points, and the 
covering defined by the rational function tp(x) is not a Belyi map. Then one of the 
following statements must hold: 

(i) Two of the three exponent differences of E\ are equal to 1/2; or 

(ii) Ei has a basis of solutions consisting of algebraic functions of z. 

Proof. Let D = degc/?, and A = {0, 1, oo} C P*. Since tp: P* -> ¥] is not a Belyi 
map, there is a branching point Pq that does not lie above A. By part (a) of 



Lemma 2.2 there are at least D + 3 distinct points above A. At most 3 of them 
can be singularities of the pulled-back equation, because Po will be a singularity by 
Lemma |2.1[ 6). Therefore there are at least D ordinary points above A. 

One or more of the 3 exponent differences of E\ must be of the form 1/fc for an 
integer fc 1, because only then ordinary points occur above A by Lemma |2.1[ c). 
Above a point of A with the exponent difference 1/fc, there may be at most D/k 
ordinary points. Let M denote the number of restricted exponent differences of E\. 
There are three possibilities: 

o M = 1. One must have k = 1, and by Lemma |2.1|[e), this point is not 



a relevant singularity for E\. Let m denote the number of distinct points 
above the two (generally) relevant singularities of E\, If m — 2, the covering 
is cyclic (i.e., Mobius-equivalent to ip(x) = x D ). If m = 3, there is only one 
branching point not above the relevant singularities of E\, by Lemma |2.2| ^fe) 
basically. Hence <p is a Belyi covering for m ^ 3. If m > 3, the pulled-back 
equation will have more than 4 singularities. 

M = 2. The exponent differences will be 1/fc, l/£ with k,£ positive integers 
and D/k + D/£ ^ D. One must have 1/fc, l/£ = 1/2, which is case (i). 
M = 3. The exponent differences will be 1/fc, l/£, 1/m with k,£, m positive 
integers and D/k + D/£ + D/m > D, i.e., l/k + l/i+l/m > 1. The subcase 
1/k + l/£ + 1/m = 1 can be ruled out, since even if the points above A 
are optimally arranged, there will be fewer than D ordinary points above 



A, contradicting Lemma 2.2 'a). It is known OHO] that in the subcase 



l/fc + l/£ + l/m>l, the equation E\ has only algebraic solutions. □ 



Remark. In case (i), the projective monodromy group of E\ is generally an infi- 
nite dihedral group. As we recall in f|5j the possible projective monodromies in 
case (ii) are: a finite cyclic, a finite dihedral, A4 (tetrahedral), S4 (octahedral) or 
(icosahedral) groups. If M = 1, the monodromy is generally an infinite cyclic 
group. Gauss-to-Heun transformations with (finite or infinite) cyclic or dihedral 
monodromy are going to be considered throughly in a separate article |32| . 
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Exponent differences Deg. Branching pattern Covering characterization, 



hyperg. Heun D above singularities composition 

a, ft, 7 a, a, 2a, 2*y 2 2 = 1+1 = 1+1 #32, Fx, F*, indecomposable 

1/2, a, P a, a, 2a, 4/3 4 [2] 2 = 4 = 2+l+l H35, Fa, Fi , 2 ■ 2 

a, 3a, P, 3P [2] 2 = 3+1 = 3+1 H 47 , F' A , F$ , indecomposable 

2a, 2a, f3, 3/3 [2] a = 3+1 = 2+2 no covering, N 27 

2a, 2a, 2/3, 2/3 [2] 2 = 2+2 = 2+2 H 3 x, F 3 , F 4 * , 2 X 2 

1/2, a, 2a, 3/3 3 [2]i+l = 2+1 = 3 H34, F 2 , F$, indecomposable 

1/3, q, /3 a, 2a, J, 2/3 ~~ 3 [3]i = 2+l = 2+l H 34 , Fj', F 3 *, indecomposable 

a, a, a, 3/3 [3]i=3=l+l+l f/33, indecomposable 

Table 1 . Possible branching patterns of hypergeometric-to-Hcvm 
transformations with 2 or 3 free parameters. 



3. Main result: Generation and classification 
Here we present the method and the results of classification of Gauss-to-Hevm 



transformations with at least one free parameter. Following part (c) of Lemma 2.1 



we restrict m € {0, 1, 2} local exponent differences of the general hypergeometric 



equation ( 1.1 ) to the reciprocals of integers k ^ 1. Thereby we have M = 3 — m free 
parameters. Basically, the free parameters are the unrestricted exponent differences. 

We ignore the cases (considered in |32j l when an exponent difference is restricted 
to 1 at a non-logarithmic singularity, or when two exponent differences are restricted 
to 1/2, as then the hypergeometric equation has cyclic or dihedral monodromy. 
Apart from this, Tables [TJ [2j [3] below give a full list of Gauss-to-Heun pull-back 
transformations with a free parameter in terms of the exponent differences (in the 
first two columns), the degree and the branching pattern of the pull-back covering 
(in the next two columns) among the entries where a covering is indicated by the H- 
notation in the last column. Table [4] gives a full list of the encountered Belyi maps 
(up to Mobius transformations and complex conjugation), and the introductory 
Figure [T] depicts the dessins d'enfant of those Belyi maps. The parallel article [30] 
identifies the pulled-back Heun equations in detail, and gives a representative list 
of transformation formulas between hypergeometric and Heun functions. 

We proceed to explain the results and notation in Tables [T] [4] Let u) denote a 
primitive cubic root of unity, say to — exp(27ri/3). In particular, u> 2 + ui + 1 = 0. 

The pull-back transformations from a hypergeometric equation E\ to a Heun 
equation E2 are classified and demonstrated in the following four steps. They 
parallel the principal steps (i)-(iii) outlined in the introduction, with the only 
difference that Step (?) is split into two steps. 

Step 1 is determination of possible restrictions on the exponent differences of E\ 
and the degree of the pull-backs. This step is elaborated in §3.1| The restrictions on 
the exponent differences determine the type of possible branching patterns, which 
is by definition an unordered list of the integers k ^ 1 that determine the restricted 
exponent differences 1/k. The following list of types is obtained: 

(3.1) (), (2), (3), (2,3), (2,4), (2,5), (2,6), (3,3), (3,4), (4,4). 

The first type () means no restrictions on the parameters of E\. We skipped 
the types (1) and (2,2); they are considered in [32] as mentioned. The types are 
indicated by the exponent differences of Ex in the first columns of Tables [Tl |3j and 
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Exponent differences Deg. Branching pattern Covering characterization, 
of the Heun equation D above singularities composition 



a, a, a, 9a 


12 [2] 6 = [3] 4 = 9+1+1+1 


H u 3 C -4 


a, a, 2a, 8a 


[2] 6 = [3] 4 = 8+2+1+1 


Hi, F23, F34, 2-2-3 


a, a, 3a, 7a 


[2] 6 = [3] 4 = 7+3+1+1 


no covering, Ni 


a, 2a, 2a, 7a 


[2] 6 = [3] 4 = 7+2+2+1 


no covering, 7V 2 


a, a, 4a, 6a 


[2] 6 = [3] 4 = 6+4+1+1 


no covering, iV 3 


a, 2a, 3a, 6a 


[2] 6 = [3] 4 = 6+3+2+1 


H 3 , F 27 , F 3 * 3 , 3-4, 4-3 


2a, 2a, 2a, 6a 


[2] 6 = [3] 4 = 6+2+2+2 


no covering, A^ 4 


a, a, 5a, 5a 


[2] 6 = [3] 4 = 5+5+1+1 


Hi, F 2i , F 3 * 2 , 2 H ■ 6 


a, 2a, 4a, 5a 


[2] 6 = [3] 4 = 5+4+2+1 


no covering, TVs 


a, 3a, 3a, 5a 


[2] 6 = [3] 4 = 5+3+3+1 


no covering, 


2a, 2a, 3a, 5a 


[2] 6 = [3] 4 = 5+3+2+2 


no covering, N7 


a, 3a, 4a, 4a 


[2] 6 = [3] 4 = 4+4+3+1 


no covering, Ns 


2a, 2a, 4a, 4a 


[2] 6 = [3] 4 = 4+4+2+2 


H 5 , F22, F£ u 2 -3c -2, 2x2-3 


2a, 3a, 3a, 4a 


[2] 6 = [3] 4 = 4+3+3+2 


no covering, Ng 


3a, 3a, 3a, 3a 


[2] 6 = [3] 4 = 3+3+3+3 


He, 3c • 4, 2h ■ 2h ■ 3c 


1/3, a, a, 8a 


10 [2] B = [3] 3 +1 = 8+1+1 


Hi, indecomposable 


1/3, a, 2a, 7a 


[2] 5 = [3] 3 +l = 7+2+1 


Hg, F21, F2S, indecomposable 


1/3, a, 3a, 6a 


[2] 5 = [3] 3 +l = 6+3+1 


no covering, iVio 


1/3, 2a, 2a, 6a 


[2] 5 = [3] 3 +l = 6+2+2 


no covering, Nn 


1/3, a, 4a, 5a 


[2] 5 = [3] 3 +l = 5+4+1 


Hg, Fig, F 2 * 9 , indecomposable 


1/3, 2a, 3a, 5a 


[2] 5 = [3] 3 +l = 5+3+2 


H10, F26, F£ , indecomposable 


1/3, 2a, 4a, 4a 


[2] 5 = [3] 3 +l = 4+4+2 


no covering, N12 


1/3, 3a, 3a, 4a 


[2] B = [3] 3 +1 = 4+3+3 


no covering, JV13 


1/2, a, a, 7a 


9 [2] 4 +l = [3] 3 = 7+1+1 


Hu, indecomposable 


1/2, a, 2a, 6a 


[2] 4 +l = [3] 3 = 6+2+1 


H12, F20, F21, 3 • 3 


1/2, a, 3a, 5a 


[2] 4 +l = [3] 3 = 5+3+1 


H13, Fis, F£e, indecomposable 


1/2, 2a, 2a, 5a 


[2] 4 +l = [3] 3 = 5+2+2 


no covering, iVi 4 


1/2, a, 4a, 4a 


[2] 4 +l = [3] 3 = 4+4+1 


no covering, A^is 


1/2, 2a, 3a, 4a 


[2] 4 +l = [3] 3 = 4+3+2 


H14, F25, F25, 3 • 3 


1/2, 3a, 3a, 3a 


[2] 4 +l = [3] 3 = 3+3+3 


no covering, Nie 


2/3, a, a, 6a 


8 [2] 4 = [3] 2 +2 = 6+1+1 


H u , F 14 , 2-4 


2/3, a, 2a, 5a 


[2] 4 = [3] 2 +2 = 5+2+1 


Hie, Fir, indecomposable 


2/3, a, 3a, 4a 


[2] 4 = [3] 2 +2 = 4+3+1 


no covering, A^i7 


2/3, 2a, 2a, 4a 


[2] 4 = [3] 2 +2 = 4+2+2 


no covering, A^is 


2/3, 2a, 3a, 3a 


[2] 4 = [3] 2 +2 = 3+3+2 


Hu, F13, 2-4 


1/3,1/3, a, 7a 


[2] 4 = [3] 2 +l+l = 7+l 


His, indecomposable 


1/3, 1/3, 2a, 6a 


[2] 4 = [3] 2 +l+l = 6+2 


H19, F w , F2* u 4 S -2, 2-4 


1/3, 1/3, 3a, 5a 


[2] 4 = [3] 2 +l+l = 5+3 


no covering, Nig 


1/3, 1/3, 4a, 4a 


[2] 4 = [3] 2 +l+l = 4+4 


H20, Fi 5 , Fi , 4-2, 2 H -4 A 


1/2,1/3, a, 6a 


7 [2] 3 +l = [3] 2 +l = 6+l 


H21, indecomposable 


1/2,1/3, 2a, 5a 


[2] 3 +l = [3] 2 +l = 5+2 


H22, Fu, Fis, indecomposable 


1/2, 1/3, 3a, 4a 


[2] 3 +l = [3] 2 +l = 4+3 


H23, F12, Fig, indecomposable 



Table 2. Possible branching patterns for pull-back transforma- 
tions from £7(1/2, 1/3, a) to a Heun equation, of degree D > 7. 
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Exponent differences 



hyperg. 



Heun 



Deg. Branching pattern 
D above singularities 



Covering characterization, 
composition 



1/2,1/3, a 



1/3,2/3, a, 5a 6 
1/3,2/3, 2a, 4a 
1/3,2/3, 3a, 3a 
1/3,1/3,1/3, 6a 
1/2,1/2, a, 5a 
1/2,1/2, 2a, 4a 
1/2,1/2, 3a, 3a 
1/2,2/3, a, 4a 5 
1/2,2/3, 2a, 3a 
1/2,1/3,1/3, 5a 
1/2,1/2,1/3, 4a 4 



3 = [3]i+2+l = 5+l 
3 = [3] i+2+l = 4+2 
3 = [3]i+2+l = 3+3 
3 = [3]i+l+l+l = 6 
2+1+1 = [3] 2 = 5+1 
2+1+1 = [3] 2 =4+2 
2+1+1 = [3] 2 = 3+3 
2 +l = [3]i+2 = 4+l 
2 +l = [3]i+2 = 3+2 
2 +l = [3]i+l+l = 5 
i+l+l = [3]i+l = 4 



#24, Fg, indecomposable 
H 25 ,F 10 , 3-2 
no covering, A/20 
#38, 3c • 2 
#26, indecomposable 
#27, F r , F? 3 , 2-3 
#28, Fq, F 12 , 2h ■ 3c 
#29, F$, indecomposable 
#30, F 5 , indecomposable 
#37, indecomposable 
#36, indecomposable 



1/2,1/4, a 



a, a, a, 5a 
a, a, 2a, 4a 
a, a, 3a, 3a 
a, 2a, 2a, 3a 
2a, 2a, 2a, 2a 
1/2, a, a, 4a 
1/2, a, 2a, 3a 
1/2, 2a, 2a, 2a 
1/4,1/4, a, 5a 
1/4, 1/4, 2a, 4a 
1/4, 1/4, 3a, 3a 
1/2,1/4, a, 4a 
1/2,1/4, 2a, 3a 
1/2,1/2, a, 3a 
1/2, 1/2, 2a, 2a 
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4 = [4] 2 = 5+1+1+1 
4 = [4] 2 = 4+2+1+1 
4 = [4] 2 = 3+3+1+1 
4 = [4] 2 = 3+2+2+1 
4 = [4] 2 = 2+2+2+2 
3 = [4]i+2 = 4+l + l 
3 = [4]!+2 = 3+2+l 
3 = [4]i+2 = 2+2+2 
3 = [4]i + l + l = 5+l 
3 = [4]! + l + l=4+2 
3 = [4]i + l + l = 3+3 
2+l = [4]i + l=4+l 
2+l = [4]x + l = 3+2 
i+l+l = [4]i= 3+1 
i+l+l = [4]i= 2+2 



no covering, A/21 

#40, Fg, F24, 2-2-2 

H20, F$, F 2 * 3 , 4-2, 2 H -U 

no covering, A/22 

#41, -F7, F 2 * 2 , 2x2x2 

no covering, N23 

#25, Fq, F* 5 , 3 • 2 

no covering, A/24 

#42, indecomposable 

no covering, A/23 

#43, -F14, Fi 4 , 2h • 3 

#44, indecomposable 

#29, .Fi'i, Fig, indecomposable 

#36, indecomposable 

#35, Fj'o, F 7 *, 2 h -2 



1/2,1/5, a 



1/5, a, a, 4a 
1/5, a, 2a, 3a 
1/5, 2a, 2a, 2a 
1/2, a, a, 3a 
1/2, a, 2a, 2a 



3 = [5] 1 +1 = 4+ 1 + 1 #42 , indecomposable 

3 = [5] l+l = 3+2+1 #24, -P15, Fie, indecomposable 

3 = [5] i+l = 2+2+2 no covering, /V 25 

2+1 = [5]i = 3+1+1 #37, indecomposable 

2+1 = [5]i = 2+2+1 #45, F{ 2 , Fn, indecomposable 



1/2,1/6, a 



a, a, a, 3a 
a, a, 2a, 2a 



3 = [6] 
3 = [6] 



1 = 3+1+1+1 
1 = 2+2+1+1 



#38, 3c • 2 

#39, F13, F17, 3 • 2, 2h ■ 3 



1/3, 1/3, a 



a, a, a, 3a 
a, a, 2a, 2a 
1/3, 1/3, a, 3a 
1/3,1/3, 2a, 2a 



2 = [3] 2 = 3+1+1+1 no covering, A/ 2 6 

2 = [3] 2 = 2+2+1+1 #28, 2 -3c 

i+l = [3] i+l = 3+1 #46, F'l , F g * , indecomposable 

i+l = [3]i + l = 2+2 #47, F3', F 8 *, indecomposable 



1/3, 1/4, a 
1/4,1/4, a 



1/3, a, a, 2a 
a, a, a, a 



= [4] 



[4]! = 2+1+1 

L = 1+1+1+1 



#36, indecomposable 
#48, 2h ■ 2 



Table 3. The other possible branching patterns of Gauss-to-Heun 
transformations with one free parameter. 
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Id Deg. Branching pattern A rational expression for ifi(x) 

#1 12 [2] 6 = [3] 4 = 9+1+1+1 64 a; 3 (a; 3 - l) 3 /(8a; 3 -9) 

#2 [2] 6 = [3] 4 = 8+2+1 + 1 27a; 2 (a; 2 -4)/4(a; 4 --4a; 2 + l) 3 

#s [2] 6 = [3] 4 = 6+3+2+1 27 (x- l) 3 (2x - 3) 2 (x + 3)/4x 3 (x 3 - 6a; + 6) 3 

# 4 [2] 6 = [3] 4 = 5+5+1 + 1 1728a: 5 (a; 2 -lla;-l)/(a; 4 -12a; 3 + 14a; 2 + 12a; + l) 3 

#5 [2] 6 = [3] 4 =4+4+2+2 27 x\x 2 - l) 2 /4 (a; 4 - x 2 + l) 3 

#6 [2] 6 = [3] 4 = 3+3+3+3 -64a; 3 (a; 3 - l) 3 /(8x 3 + l) 3 

H 7 10 [2] 5 = [3] 3 +l = 8+1+1 -4(x + 2)(a; 3 + 3x + 2) 3 /27(3a; 2 -2a; + ll) 

H 8 [2] 5 = [3] 3 +l = 7+2+1 4(a; + 4)(a; 3 -6a;-2)727(3a; + 4) 2 (4a;-ll) 

#9 [2] 8 = [3] 3 + l = 5+4+1 -(8a; - l)(8a; 3 + 87a; 2 + 96a; - 64) 3 /2 3 3 12 a; 4 (a; + 10) 

H w [2] 6 = [3] 3 + l = 5+3+2 -(a; - 3) (81a; 3 - 9a,- 2 - 53a; - 27) 3 /2 14 3 a; 3 (9a; + 5) 2 

#n 9 [2] 4 +l = [3] 3 = 7+1 + 1 4(a; 3 + 4x 2 + 10a; + 6) 3 /27(4a; 2 + 13a; + 32) 

#12 [2] 4 +l = [3] 3 = 6+2+1 27a,- 2 (a;-3)/4(x 3 -3a; 2 + l) 3 

H 13 [2] 4 +l = [3] 3 = 5+3+1 -25(5x 3 + 45a; 2 + 39a; - 25) 3 /2 14 3 3 a- 3 (3a; + 25) 

#14 [2] 4 +l = [3] 3 =4+3+2 27 a; 3 (3a; - 4) 2 /4 (a; 3 -3a; -4) 3 

H 15 8 [2] 4 = [3] 2 +2 = 6+1+1 64a; 2 (a; 2 - l) 3 /(8a; 2 -9) 

H le [2] 4 = [3] 2 +2 = 5+2+1 -4a; 2 (a; 2 + 8a; + 10) 3 /27(2a; + l) 2 (4a; + 27) 

H 17 [2] 4 = [3] 2 +2 = 3+3+2 -64a; 2 (a; 2 - l) 3 /(8x 3 + l) 3 

#18 [2] 4 = [3] 2 + l+l = 7+l -(a; 2 -13a; + 49)(a; 2 -5x+l) 3 /2 6 3 3 a; 

H w [2] 4 = [3] 2 + l+l = 6+2 -64a; 2 /(a; 2 - l) 3 (a; 2 -9) 

H 20 [2] 4 = [3] 2 + l+l = 4+4 16a; 3 (2a; + l)(a;-4)/(a; 2 -2a;-2) 4 

# 2 i 7 [2] 3 + l = [3] 2 + l = 6+l 4(a;-l)((l + 2cj)a; 2 -3a;-a;) 3 /(4-(l + 3a;)a;) 

# 22 [2] 3 +l = [3] 2 + l = 5+2 a;(2a; 2 - 35a; + 140) 3 /108(14a; - 125) 2 

#23 [2] 3 + l = [3] 2 + l = 4+3 -(a; + 27)(16a; 2 + 80a; - 243) 3 /2 2 3 3 7 7 a; 3 

# 24 6 [2] 3 = [3]i+2+l = 5+l -a; 2 (a;-3)(a; + 5) 3 /64(3x + 16) 

# 25 [2] 3 = [3]i+2+l = 4+2 -4x 3 (x- l) 2 (a; + 2)/(3a;-2) 2 

#26 [2] 2 +l+l = [3] 2 = 5+l 1728a;/(a; 2 + 10a; + 5) 3 

#27 [2] 2 + l+l = [3] 2 =4+2 27a; 2 /4(a; 2 - l) 3 

#28 [2] 2 + l+l = [3] 2 = 3+3 (a; 2 + 6a; - 3) 3 /(a; 2 - 6a; - 3) 3 

#29 5 2+2+1 = 3+2 = 4+1 4a; 2 (a; + 5) 3 /27(5x + 27) 

# 30 2 +2+1 = 3+2 = 3+2 a; 2 (9a;-5) 3 /4(5a;-l) 2 

# 3 i 4 2+2 = 2+2 = 2+2 -4x 2 /(a; 2 - l) 2 

# 32 2 2 = 2= 1+1 a; 2 

# 33 3 3 = 3= 1+1+1 a; 3 

# 34 3 = 2+1 = 2+1 a;(4a;-3) 2 

H 35 4 4 = 2+2 = 2+1+1 4x 2 (l~x 2 ) 

# 36 4 = 3+1 = 2+1+1 -x 3 (3a; + 4) 

# 37 5 5 = 2+2+1 = 3+1+1 (a; 2 + llx + 64)(x + 3) 3 /2 6 3 3 

# 38 6 [2] 3 =6 = 3+1+1+1 4a; 3 (1 - a; 3 ) 

#39 6 [2] 3 =6 = 2+2+1+1 a; 2 (4a; 2 - 3) 2 

# 40 8 [2] 4 = [4] 2 =4+2+1 + 1 4a; 2 (a; 2 - 2)/(x 2 - l) 4 

#4i [2] 4 = [4] 2 = 2+2+2+2 -4x 4 /(x 4 ~l) 2 

#42 6 [2] 3 = 4+1+1 = 5+1 -(a;-l) 4 (a; 2 -6a; + 25)/256a; 

# 43 [2] 3 = 4+1+1 = 3+3 27a;(a;-4)/4(a; 2 -4a; + l) 3 

#44 5 4+1=4+1 = 2+2+1 a;(a;-l-2i) 4 /((l + 2i)a;-l) 4 

#45 5 = 2+2+1 = 2+2+1 x \x 2 - 5a; + 5) 2 /4 

#46 4 3+1 = 3+1 = 3+1 a; 3 (a- + 2)/(2a; + l) 

#47 3+1 = 3+1 = 2+2 64a;(a;-l) 3 /(8x-9) 

H i8 4 = 4=1+1+1+1 a; 4 



Table 4. The Belyi coverings appearing in Gauss-to-Heun pull- 
backs, up to Mobius transformations. 
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the whole Table [2] is devoted to the type (2,3). The entries of different types are 
separated by horizontal lines. The pull-back degree is given in the third columns 
of Tables [T] [3j and the second column of Table [2] The maximal degree is 12. It 
occurs for the type (2,3) only. 

Step 2 is determination of possible branching patterns. The method is explained 
in §3.2| The result is presented by the fourth columns of Tables [T] [3j and the third 
column of Table [2] Generally, we indicate a branching pattern by an (unordered) 
list of three unordered partitions of its degree D, separated by the equality signs. 
The partitions specify the branching indices in each of the three branching fibers 
of a Bclyi covering. Besides, we use the abbreviation [k) n for a partition block 
k + ■ ■ ■ + k (n times). In Tables [I] [2] [3j the symbol [k] n specifically means presence 
of n points of E 2 with the branching order k above a singular point of Ei with the 
the exponent difference 1/k. By part (c) of Lemma |2.1| each of the n points will 
be either ordinary or an irrelevant singularity for E 2 . By the described convention, 
the branching patterns for pull-back transformations with M free parameters have 
3 — M numbers (i.e., branching orders) inclosed in square brackets, and exactly 4 
non-bracketed numbers representing the 4 singular points of E 2 ■ 

In total, we get a list of 89 branching patterns, though some of the patterns differ 
only by the square-brackets specification of ordinary points of E 2 . For example, 
two degree 3 branching patterns in Table [l] are the same, leading to the same cubic 
covering H34 (identified in the last column). The exponent differences of E 2 are 
determined by E\ and the branching pattern, and are given by the second columns 
of Tables [TJ [3] and the first column of Table [2] 

Step 3 is computation of the Belyi coverings ip : P* — > . Generally, computa- 
tion of Belyi maps with a given branching pattern is a difficult problem. However, 
the maximal degree implied by the possible branching patterns is just 12. With the 
aid of modern computer algebra systems this problem is tractable for coverings of 
degree 12 or less, even using a straightforward Ansatz method with undetermined 
coefficients. Most of the Belyi maps are actually known in the literature, if only 
because the Belyi maps of the type (2,3) occur in Herfurtner's list [TU] of elliptic sur- 
faces with four singular fibers. Specifically, the J{X, y)-expressions in JTOJ Table 3] 
are homogeneous expressions of the Belyi maps H\ , . . . , up to Mobius trans- 
formations. Moreover, the same coverings basically appear in pull-backs between 
hypergeometric equations, because a free parameter can always be specialized so 
to reduce the Heun equation E 2 to a hypergeometric (or simpler) equation. 

The computational issues of Step 3 are discussed in |3.2| Complementarily, 
f|5] presents an elegant approach to show non-existence of Belyi maps with many 
branching patterns. The full list of computed Belyi maps is given in Table |4j and 
further commented in Sj4] The last columns of Tables [T] [2] [3] identify the Belyi 
map for each possible pull-back transformation. These columns also specify the 
Coxeter decompositions [7 and divisible tilings [5] for the Schwarz maps associ- 
ated to the pulled-back Heun's equation E 2 (by various -F-numbers, as explained 
in £ 4.3 ) , and describe composite transformations by product expressions indicating 
degrees of occurring indecomposable transformations. The product notation has to 
be followed from right to left to trace the composition from the starting hypergeo- 
metric equation. The factor 2h denotes quadratic Heun-to-Heun transformation 
(2.7). Here is the meaning of other indexed degrees: 3c denotes the cyclic cover- 
ing H33 with the branching pattern 3 = 3 = 1+1+1, while A a and 4b stand for the 
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coverings H 36 (4 = 3+1 = 2+1+1) and H 46 (3+1 = 3+1 = 3+1), respectively. The 
unindexed numbers 3 and 4 denote the frequent coverings H 3 4 (3 = 2+1 = 2+1) and 
i?47 (3+1 = 3+1 = 2+2), respectively. In any composition, there is exactly one fac- 
tor representing an indecomposable Gauss-to-Heun transformation; it is the first 
one from the left which is not 2#. The other factors to the right represent pull-backs 
between hypergeometric equations. The notation 2x2 indicates a composition of 
quadratic transformations that can be realized in multiple ways, possibly including 
2h] see (4.3) below for the most typical example. The compositions are considered 
more thoroughly in [4.4 and in [30, Appendix B]. 

There are 27 different branching patterns for which there is no Belyi map. The 
non-existence in all these cases can be elegantly shown by considering implied (but 
not possible) pull-back transformations between Fuchsian equations, as explained 
in <|5] The indexed TV-notation refers to Table [5] below. For each branching pattern 
except two leading to if 21 and -f/44, there is at most one covering (and one pull-back) 
up to Mobius transformations. The coverings H21 and H44 are defined, respectively, 
over Q(w) and Q(i). In either of these cases, we actually have a complex-conjugated 
pair of Belyi coverings. Table [| lists 48 different coverings, though H 2 \ and -H44 
should be properly counted twice. It is instructive to compare the branching pattern 
and the orders of vertices and cells of the dessins d 'enfant in Figure [T] In total, we 
count 61 parametric pull-backs among the entries of Tables [lj [2j [3] Of them, 28 are 
composite. Evidently, some of the 48 coverings appear in more than one pull-back. 
Accordingly, the symbol [k] n in Table [4] merely indicates presence of n points of 
branching order k in the same fiber. The coverings H 20 , H 2i , H 25 , H 2S , H 29 , -H34, 
i?35, i?37, -^38, -^42, #47 appear twice in Tables [T] [2j [3j while H%q three times. 

Step 4 is derivation of identities between standard iFi{z) and Hn(i) solutions 
of the related hypergeometric and Heun equations, with z = fix). This gives Heun- 
to-hypergeometric reduction formulas, expressing found Heun functions in terms of 
the better understood Gauss hypergeometric functions. This final step is compre- 
hensively considered in the parallel paper [30) by the same authors. In particular, 
[3UI §3] explains the technical issue of choosing the gauge prefactor 0(x) in pull-back 
transformations (1.4). The transformations without a prefactor (i.e., 6{x) — 1) are 
classified by Maier in [16]. The branching patterns for these pull-backs typically 
have a fiber with just one point, and that point is a singularity for E 2 . There are 
7 of these pull-back transformations. Their type is (), (2), (3) or (2,3), and the 
coverings are numbered consequently from H 32 to H^g. Formulas without a prefac- 
tor arise from the transformations of Tables [l] [2] realized by these coverings, except 
for the type (3) transformation with the covering H34. The well-known quadratic 
transformation (2.8) is described at the beginning of this article. 

Hereby we complete the description of four classification steps. At the same 
time, we explained the results and notation in Tablcs[I]-[4] The next two subsections 
give a methodological proof of Steps 1 and 2. Section [4] discusses computational 
issues of Step 3, composite coverings, and relations of the recorded transformations 
to the Herfurtner's list [lOj of elliptic surfaces and Fclikson's list [7] of Coxeter 
decompositions. Section [5] describes the elegant approach of proving non-existence 
of Belyi coverings with certain branching patterns, and applies it to Tables [Tf|3] and 
the Miranda-Persson list [TS] of degree 24 branching patterns. 
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3.1. Step 1: Possible restricted exponent differences and degree. We are 

looking for the Belyi coverings ip : P* —> V\ that pull-back a hypergeometric equa- 
tion E\ to Heun's equation E 2 . We assume that Ei is not specifically of the form 
E(l,a,j3) or E(l/2, 1/2, a), because then either it has a logarithmic singularity 
(if j3 7^ ±a by [29, Lemma 5.1]; it would not contribute extra non-singular points 
above {0,1, oo} C V\ necessary for new cases of Gauss-to-Heun pull-backs), or it 
has cyclic (if (3 = \a\) or dihedral monodromy as explored in |32j . 

We restrict m E {0,1,2} exponent differences of the general hypergeometric 
equation ( |1.1[ ) to t he re ciprocals of integers k > 1, and look for particular cases when 



part (c) of Lemma 2.1 allows enough non-singular points above {0, 1, oo} C P*. The 
degree of <p is denoted by D. 

First, assume that m = 0. This puts no restriction on the exponent differences 
of Ei, so generally all points above z = 0, 1, oo will be singularities of E 2 - There 
will be exactly D + 2 singular points by Lemma |2.2| and we wish the transformed 
equation to be Heun's. Hence D+2 ^ 4, so that D ^ 2. If D = 1 then <p is a Mobius 
transformation and does not alter the number of singular points, hence Ei will have 
only three. If D = 2, then ip is a quadratic covering with the branching pattern 
2 = 2 = 1 + 1. The pull-back is then the well-known quadratic transformation 



(2.8), applicable to any hypergeometric equation. We do not need to consider 
quadratic transformations (D — 2) subsequently, nor the case m = in more detail 
in the other steps. 

For m € {1, 2}, the number of non-singular points above the restricted singular- 
ities of Ei must be at least (D + 2) — 4 = D — 2. 

If m — 1, we allow two free parameters. We restrict just one exponent difference 
of Ei to equal 1/fc, with integer k > 1. The pulled-back equ ation E 2 will have 
at most \_D/k\ ordinary points above {0, 1, oo} C by Lemma 2.1 and one must 
have 



(3.2) 



> D-2. 



This leads to the Diophantinc inequality 

(3.3) £ + J>L 

For k > 1 and D > 2, we have the following possibilities: 

(3.4) (k,D)e{(2,3), (2,4), (3,3)}. 

The resulting branching patterns are of types (2), (3), according to k. 

If m = 2, we allow one free parameter. Suppose that the restricted exponent 
differences of Ei equal 1/fc, l/£, where k,£ are integers. We assume 1 < k ^ I ^ D 
without loss of generality; the last inequality allows actual utilization of the restric- 
tion l/£. The transformed equation has at most \D/k\ + [D/£\ ordinary points 
above {0, 1, oo} C P*. Similarly to the above, one must have 



(3.5) 



D 




D 


_~fc_ 


+ 


J _ 



> D-2, 



which leads to the weaker Diophantine inequality 
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Discarding k = £ = 2, we get the following possibilities for k,£ and for the upper 

bound -D max on the degree D: 

(3.7) 

(k,£,D max )e {(2,3,12), (2,4,8), (2,5,6), (2,6,6), (3,3,6), (3,4,4), (4,4,4)}. 

The resulting branching patterns are the seven types (2,3), (2,4), (2,5), (2,6), 
(3,3), (3,4), (4,4), respectively. 

3.2. Step 2: Possible branching patterns. Here we look at each possible type 
and degree, and determine all branching patterns fitting them. The constraint that 
there must be at least D— 2 ordinary points above {0, 1, oo} C P* will usually require 
taking the number of ordinary points above the points with restricted exponent 
differences is maximal, i.e., equal to \_D/k\ or \D/t\. 

The a priori possible branching patterns for the case m = 1 are straightforward 
to determine. They are listed in the fourth column of Table [l] That table is 
comparable to Table lb]. 

In the case m = 2, we start with the coverings of the type (2, 3) of the maximal 
degree D = 12, as in Table[2] There must be 12 — 2 = 10 ordinary points above the 
two singular points of £(1/2,1/3, a) with exponent differences 1/2 or 1/3; all IE- 
points in these two fibers must be ordinary, as |^12 /2J + [12/3J = 10. The third fiber 
is a partition of 12 with 4 parts. There are 15 such partitions, and they are all listed 
in the third column of Table [2j Next, there are no transformations of degree 11, 
because 11 — 2 > |11/2J + UV^J an d there would not be enough ordinary points 
in the two fibers. In a similar way, the pull-back coverings of degree D = 10, 9, 8 
or 7 must have the maximal number of ordinary points in the two restricted fibers; 
and all branching patterns consistent with this constraint are listed. The branching 
patterns of type (2, 3) continue in Table [3j The degrees D = 6, 4 require less than 
L-D/2J + [-D/3J ordinary points in the restricted fibers, and there is some choice of 
how to split a bracketed number [2] or [3] into a pair of non-bracketed numbers, 
though at least one bracketed number must remain in the two restricted fibers. For 
D = 5, there is a choice of splitting (or not splitting) the number 2 in the [3] fiber. 
In total, we get 53 branching patterns of the type (2, 3), all different. 

The other types (2, 4), (2, 5), (2, 6), (3, 3), (3, 4), (4, 4) similarly give less numer- 
ous sets of branching patterns, some of them coinciding mutually or with previously 
encountered ones. 



4. The Belyi coverings 

First, this section briefly explains computation of Belyi maps and utilizing spe- 
cialization of parametric Gauss-to-Heun transformations to transformations be- 
tween hypergeometric equations. In § §4.2|p~3"l we explain how the if-numbering of 
Table [4] comes partly from an algebraic-geometric classification of Herfurtner |10j , 
and clarify the various F-numbers in the last columns of Tables [T}|3] as representing 



Coxeter decompositions of Felikson [7] and divisible tilings of [5]. Lastly, in S4.4 
examine the composite transformations among our results. 

4.1. Computational issues. To compute the Belyi maps cp: ¥]. — > with a 
given branching pattern means to find all rational functions tp(x) such that the 
numerators of <p(x), 1 — <p(x) and the denominator of (p(x) factor according to the 
branching pattern. A straightforward Ansatz method with undetermined coeffi- 
cients can be used for low degree coverings. Modern computer algebra systems 
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(such as Maple and Mathematica) can handle the resulting systems of algebraic 
equations easily if the degree of ip(x) is 12 or less. More cannily, one can consider 
factorization of the numerators of the logarithmic derivatives of y>(x) and ip(x) — 1 
as in [261 §3]. For example, to determine Hi, one is looking for a constant c and 
monic polynomials P, Q, R of degree 4, 3, 6, respectively, such that ip(x) = cP 3 /Q 
and (p(x) — 1 = cR 2 /Q. To find these polynomials, one considers 

(4 . (p'{x) _ 3P' Q' j_ 9R {tp(x) - 1)' _ 2R' Q' j_ 9P 2 

( ' ) ip(x) ~ ~P ~ ~Q ~ ~PQ' (p(x) - 1 ~~ ~R ~ ~Q ~ R~Q' 

Zeroes of the derivatives are the branching points other than in the denomina- 
tors, and the factor 9 is determined by local consideration at x = oo. The whole 
polynomial R can be eliminated symbolically using the first identification, and the 
resulting equation system for the undetermined coefficients of P, Q is rather trans- 
parent. In general, a covering with a given branching pattern may not exist, or 
there may be several Belyi maps (up to Mobius equivalence) or even several Q/Q- 
Galois orbits of Belyi maps with the same branching pattern. The Galois action on 
the Belyi maps and their dessins denfant is of primary interest to Grothendieck's 
theory [22]. p4]. 

Less demandingly, one may notice that the free parameter of our Gauss-to-Heun 
transformations can be specialized so that to the pulled-back Heun equation has ac- 
tually less than 4 singular points. Therefore, the Belyi coverings must appear in the 
classification [29] of Gauss-to-Gauss transformations in principle, though there are a 
few infinite families of those transformations (for degenerate, dihedral, algebraic or 
elliptic hypergeometric functions). In particular, each branching pattern of Table [l] 
can be found in [29, Table 1], except for 2+2 = 2+2 = 2+2 which corresponds to the 

transformation £(1/2, 1/2, a) <— £(l,2a,2a) briefly mentioned in [2FJ P- 161]. 
The branching patterns of Tables [2] and [3] (with m = 1 free parameter) can be 
handled similarly, yielding reductions of one-parameter Gauss-to-Heun transforma- 
tions to zero-parameter pull-backs between hypergeometric functions. For example, 
the covering H27 implies the hypergeometric transformations £(1/2, 1/3, 1/2) <— 
£(1/2,1/2,2) and £(1/2,1/3,1/4) ^- £(1/2,1/2,1/2). These specializations re- 
ductions are possible whenever there is a branching point with a free exponent differ- 
ence. Among the relevant branching patterns, only the last one (4 = 4= 1+1+1+1) 
in Table [3] does not satisfy this condition. But even it represents a nominally hyper- 
geometric transformation, namely E(a, a, 1) <A- £(4a,4a, 1). Section^ gives more 
details for obtaining the list of Gauss-to-Heun pull-backs from the classification in 



In particular, the non-unique coverings £21 and £44 come from Lemma 5.3 



4.2. The Herfurtner classification. Pull-back transformations from hypergeo- 
metric equations of the form £(1/2, 1/3, a) to Heun equations have a close relation 
to elliptic surfaces over C(x) with 4 singular fibers [lOj [19]. The Belyi coverings 
z = ip(x) that induce these transformations appear as j-invariants of the elliptic 
surfaces, with z equal to J := j/1728, the traditional Klein j-invariant. 

The elliptic surfaces with 4 singular fibers are classified by Herfurtner [TO]. His 
article lists 50 configurations of singular fibers which give such elliptic surfaces, and 
for each configuration, supplies a formula J = J(X, Y) which is a projectivized 
version of z = <p(x), up to a Mobius transformation of x and a permutation of 
z = 0, 1, 00. Heun equations arise from 38 of his 50 cases, as Movasati and Reiter [T5] 
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recently observed. We adopt the enumeration of [HI Table 1], and denote these 
38 Belyi coverings of Herfurtner, which were not originally numbered, by Hi to 
i?38. The ordering is by degree in two ranges, as evident in Table [4j decreasing in 
the range Hi, . . . , H31, and increasing in the range H32, ■ ■ ■ , H^. 

By examining Tableland the upper part of Table [3j one finds that the coverings 
Hi, ... , H30 and -H36, • ■ • , H^s induce Gauss-to-Heun pull-backs of the type (2, 3) 
with one free parameter. These transformations use each of these 34 coverings ex- 
actly once, and no other coverings appear. The ordering by decreasing degree make 
the H- numbers appear ordered in Table [2j and almost ordered in the upper part of 
Table [3} By examining Table [IJ one finds Herfurtner 's coverings H31, . . . , H 35 (with 
H34 appearing twice) and a "new" covering H47. The covering H47 cannot pull- 
back £7(1/2, 1/3, a) to a Fuchsian equation with exactly 4 singularities. The pattern 
[3]i = 2+1 = 2+1 for H 3 4 cannot be refined to such a pull-back from E(l/2, 1/3, a) 
either, but this is possible for the other H34 parsing [2]i+l = 2+1 = 3. This explains 
why H34 appears in Herfurtner's list once. 

Some of Herfurtner's coverings additionally induce one-parameter Gauss-to-Heun 
transformations of types (2,4), (2,5), etc., as evident in Table [3j But 10 extra 
coverings appear in the latter sections of that table; they have no interpretation 
in terms of elliptic surfaces. We denote them H$g, . . . , i/48, ordered somewhat 
arbitrarily in the lower part of Table [4] The covering i/47 induces transformations 
of the types (2) and (3, 3). 

4.3. Coxeter decompositions. Recall that a Schwartz map for an second order 
differential equation in the complex domain is a map C — > C defined as the ratio of 
a pair of independent solutions of the differential equation [3] . Consider a hypergeo- 
metric equation with real exponent differences (a, /3, 7) satisfying ^ a, /3, 7 < 1. 
The image of the upper half plane under its Schwarz map is a curvilinear Schwarz 
triangle; the sides are line or circle segments, and the angles are equal to ira, tt(3, 777. 
Similarly, consider a Heun equation with real exponent differences (a,/3,j,5) sat- 
isfying ^ a, ,8,7, 6 < 1. The image of the upper half plane under its Schwarz 
map is a curvilinear Schwarz quadrangle, with the same kind of sides, and angles 
are equal to 77a, it/3, 777, ttS. 

It was noticed by Hodgkinson [TTJ [T^] that if the covering tp(x) of a pull-back 
transformation between hypergeometric equations is defined over R, the analytic 
continuations of their solutions according to the Schwarz reflection principle are 
compatible. In consequence, the covering tp (of degree D, say) will induce a sub- 
division of a Schwarz triangle of the pulled-back hypergeometric equation into 
D Schwarz triangles of the original hypergeometric equation. Examples of such 
subdivisions are given in [551 Figure 1]. 

Similarly, suppose we have a Gauss-to-Heun transformation defined over R. In 
particular, the fourth singular point x = t is real. Then the analytic continua- 
tions of the hypergeometric and Heun solutions according to the Schwarz reflection 
principle are compatible, and the covering ip (of degree D) will induce a subdivi- 
sion of a Schwarz quadrangle of the Heun equation into D Schwarz triangles of the 
hypergeometric equation. 

In the context of hyperbolic geometry, the possible subdivisions of curvilinear 
quadrangles (or triangles) into curvilinear triangles have been classified by Felik- 
son [7]; they are called Coxeter decompositions. The triangles have angles rra, it {3, 777 
satisfying a + /3 + 7 < 1. The Coxeter decompositions with a free (angle) parameter 
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are depicted in Figures 10, 11, 14 in 7J. The subdivisions of Schwarz quadran- 
gles into Schwarz triangles induced by our Gauss-to-Heun transformations defined 
over K have the same shape. In Tables [T}|3j 

o the notation Fk refers to the fcth subdivision picture in [7J Figure 14] ; these 
subdivisions are applicable to Gauss-to-Heun pull-backs of the type (2,3); 

o F' k similarly refers to (7J Figure 11]; these subdivisions are applicable the 
pull-backs of the types (2), (2,4), (2,5), (2,6); 

o Fy. similarly refers to [7] Figure 10]; these subdivisions are applicable the 
pull-backs of the type (3) or (3, 3). 

Figure a) depicts the Coxeter dexomposition F[ 3 of a quadrangle with the 
angles ira, na, 2na, 2ira into 6 hyperbolic triangles with the angles 7r/2, 7r/6, ira. It 
gives a decomposition of a Schwarz quadrangle for HE(a, a, 2a, 2a) into Schwarz 
triangles for £'(1/2, 1/6, a) induced by the type (2,6) transformation with the cov- 
ering H 39 . The Schwarz reflection principle is applied to a few edges intersecting 
at a common vertex. The decompositions 3 • 2 and 2h ■ 3 are clearly visible in 
the Coxeter decomposition. Consequently, the picture also illustrates the decom- 
position F 2 of the same quadrangle into 3 triangles with the angles 7r/3, ira, ira, 
and the decomposition F2 of a quadrangle with the angles 7r/2, 7r/2, a, 2a. Both 
decompositions are induced by the cubic covering H34. The factor 2# represents a 
Schwarz reflection between two smaller quadrangles. 

Figure ^ b) is not a quadrangle, of course. But it contains two Coxeter decom- 
positions for Gauss-to-Heun transformations of the type (3,3). If we remove the 
upper black triangle, we get the decomposition F 3 of a quadrangle with the angles 
7r/3, 7r/3, 2na, 2ira. If the left white triangle is removed, the decomposition F'l of 
a quadrangle with the angles 7r/3, 7r/3, ira, 37ra is obtained. The coverings are #47 
and H^e, respectively. 

Similarly, Figure[2]jc) includes all Coxeter decompositions for the Gauss-to-Heun 
transformations of the types (2,4) and (2,5). Here we identify the quadrangles 
(and the corresponding Belyi coverings) for the Coxeter decompositions Fg to F[ 2 , 
respectively: 

ABCF(H 25 ), ABFH{H iX ), ABDF (H 20 ), BDFH (H i0 ), 
ABML (H 35 ), ABCL (H 29 ), OCEG (#«). 

The quadrangles (and coverings) for the Coxeter decompositions F[ 4 and F[ 5 are 
KCFH (#43) and OGEE (H 2i ), respectively. 

Finally, Figure [2^ c) includes all Coxeter decompositions for the Gauss-to-Heun 
transformations specifically of the type (2,3). They are numbered from F5 to F 27 
by Figure 10]. Here are their quadrangles (and coverings), respectively: 



AOEX (H 30 ), 
AOEP (H 25 ), 
APER (H 20 ), 
AXEB {H 12 ), 



AXEZ (#»,), 
AXEQ {H 22 ), 
APEQ (H 19 ), 
APEB (H a ), 



AXEY(H 27 ), 
AXER (H 23 ), 
AOEG (H 16 ), 
ACEF (H 5 ), 
AXEC (H u ), 



AFOY(H 29 ), 
AFEO (H 17 ), 
AXED (H 13 ), 
ABEG (H 2 ) , 
APEC (H 10 ), 



AFOQ (H 24 ), 
AFOB (H 15 ), 
APED (H 9 ), 
ADEG {H^, 
ACEG (H 3 ). 



In total, there are (27 - 4) + (15 - 5) + (4 - 2) = 35 subdivisions F k , F' k , F' k ' 
representing Gauss-to-Heun transformations with exactly one parameter. 
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(6) A 



id) 



Figure 2. Coxeter decompositions for the parametric Gauss-to- 
Hcun transformations defined over R 

The subdivisions for the Gauss-to-Heun transformations with 2 or 3 parameters 
are the following: 



o The Coxeter decomposition for quadratic transformation (2.8) is repre- 
sented by a single Schwarz reflection. It can be discerned in many places in 
Figure [2] for example as the quadrangle OYCZ in picture (d). It appears 
several times in Felikson's figures, in particular as Fi=F[=F[' . 

o There are two degree 3 decompositions F2=F' 2 and F 2 . They are both 
represented by the covering ^34, as we mentioned discussing picture (a). 
The other cubic transformation (with the covering ^33) is not defined over 
K in the normalized form [SUl §4-4.4] but over Q(w), hence there is no 
Coxeter decomposition for it. 



There are three degree 4 decompositions, F3 



F' A and F 4 =i^. They 



can be discerned, for example, as the following quadrangles in picture (c), 
respectively: OBCF (H 31 ), OABC (H i7 ), OCEF (H 35 ). 
Whether a Gauss-to-Heun transformation is realized by a Coxeter decomposi- 
tion, is determined by a close inspection in Step 4 of Sj3] A necessary and sufficient 
condition is that the Belyi covering has to be defined over M after a normaliza- 
tion (by Mobius transformations) that locates 3 of the 4 singular points of Heuns 
equation as x = 0, x = 1, x = 00. In particular, the fourth singular point x = t 
has to be real, though this is not a sufficient condition. For example, a proper 
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normalization of for the type (4,4) transformation is 8ix(x 2 — l)/(x + i) 4 . 
This gives t — — 1, but the covering is not defined over HL There is one other 
example of this type: a proper normalization of i?28 f° r a type (3, 3) pull-back 
is 3(1 + 2uj)x 2 (x 2 — l)/(x 2 + lj) 3 . On the other hand, a proper normalization of 
the same H28 for a type (2,3) pull-back is defined over Q(v3), giving the Cox- 
eter decomposition Fq. There are two different Coxeter decompositions for each of 
the following coverings: i?20, ^24, H25, ^29, H34, H^, H^. Comparison of our 
classification and Felikson's list [7] provides a useful mutual confirmation. 

The considered Coxeter decompositions are parametrized, in that one or more 
of the triangular vertex angles are free to vary. For somewhat larger real values 
of the free parameter(s), the Coxeter decompositions are transfigured to spherical 
geometry of the Riemann sphere (if angles larger than tt are allowed) , as subdivisions 
of spherical quadrangles into spherical triangles with the angles satisfying an + 
/?7r + 777 > 7r. Most of the Coxeter decompositions can be transfigured to the plain 
Euclidean geometry (where a7r + /?7r + = %) as well. The exceptions are -F14, 
-F167 ^20; -F27, Fg, F4, for which the quadrangles degenerates to flat triangles. 

Broughton et al. 5 classify similar geometric objects: divisible tilings of the 
hyperbolic plane. Compared with Felikson's pictures, divisible tilings form a proper 
subset of Coxeter decompositions. The condition for a Coxeter decomposition to be 
a divisible tiling is that the quadrangle angles be equal to rr/fc, with k an integer. 
In general Coxeter decompositions, rational multiples of tt are also allowed. The 
one-parameter divisible tilings relevant here are depicted in [SJ Table 6.6]. There 
arc 34 of them; the first 6 correspond to Gauss-to-Heun transformations with 2 
or 3 parameters. Divisible tilings are indicated in Tables [T|-|3] by the notation 
Fj , . . . , i 7 ^, where the subscripts refer to the numbering in [SJ Table 6.6]. There 
are 35 — (34 — 6) = 7 relevant Coxeter decompositions with one parameter that are 
not divisible tilings; they all have the angle 27r/3. 

4.4. Composite transformations. The composite Gauss-to-Heun transforma- 
tions can be inductively deduced from a smaller set of pull-back transformations 
among hypergeometric and Heun functions. Due to the associativity of the com- 
position operation, one can always decompose a Gauss-to-Heun transformation as 
a product of the following: 

o A possibly composite Gauss-to-Gauss transformation with a free parame- 
ter, excluding Mobius transformations and from E(l, a, a) or E(l/2, 1/2, a) 
for the purposes of this article. This could be the quadratic transforma- 



tion (2.6) and one of 6 classical transformations (of degrees 3, 4 and 6) 
worked out by Goursat [5] and listed in [2HJ Table 1]. 
o An indecomposable Gauss-to-Heun transformation with at least one free 



parameter. This could be the quadratic transformation (2.8 1; one of 4 
other indecomposable transformations of Table [l] or an indecomposable 
transformation of Table [2] or [3] of degree at most 6, possibly fitting a Gauss- 
to-Gauss or a Heun-to-Heun transformation. 

A Heun-to-Heun transformation with at least one free parameter. This 



could be the quadratic transformation (2.7), or the degree 4 composite 
transformation 

(4.2) HEQ./2, 1/2, 1/2, a) ^ HE(l/2, 1/2, a, a) BE(a, a, a, a), 
realized by the covering if 31 . See [301 §4-3] for an overview. 
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Figure[3]graphically depicts all possible compositions of the preceding three types 
(Gauss-to-Gauss, etc.) The two longest boxes, centrally placed, represent quadratic 
transformations (based on the double covering H^)- The following objects and 
information are included in the figure. 

There are 7 boxes with double edges on the left and the right sides, representing 
the classical Gauss-to-Gauss transformations. The quadratic transformation ap- 
pears as the long box in the lower part; two indecomposable transformations (of 
degree 3 or 4) appear as boxes in the central part; and the remaining four classical 
transformations (of degrees 3, 4 and 6) are represented in the upper part. Of the 
latter, only the cubic transformation is indecomposable. The transformation ap- 
pearing near the upper-right corner can be decomposed in two different ways; its 
covering does not occur in Tables [IJj3j hence it is not identified by an H number. 
These 7 boxes will be called E —> E boxes. 

The 10 other boxes represent indecomposable Gauss-to-Heun transformations. 
The quadratic transformation (2.6) is represented by the long box in the upper 
part; the four indecomposable transformations of Table [l] appear in the central 
part. The three isolated boxes near the lower right corner represent the indecom- 
posable transformations of Table [3j to each of which the quadratic Heun-to-Heun 
transformation ( |2.7[ ) can be applied. The other two lowest boxes represent trans- 
formations in Table [3] that can be composed with a specialization of the quadratic 
E — » E transformation. These 10 boxes will be called E — > HE boxes. 

The vertical lines connect E — > E and E — > HE boxes whose transformations 
can be composed (perhaps after a specialization of parameters). The composed 
coverings are labeled by H numbers on the left side of each vertical line. Relevant 
specializations of the quadratic E — > E transformation are given as well. Note that 
the specializations V = \ and q = | of the quadratic E — > E transformation are 
not given, because (as stated above) the dihedral family is not considered here. 
The number of possible compositions between an E — > E box and an E — > HE 
box depends on the number of ways to identify (without degeneracy) the exponent 
differences of the intermediate hypergeometirc equation. It is instructive to compare 
compositions of the quadratic E — > HE transformation with the two hypergeometric 
transformations coming from the coverings H47 and H34. Compositions of the 
quadratic E — > E and E — > HE transformations occur as the composite quartic 
coverings H 35 , iJ 31 in Table [T] 

The =>■ symbols outside the boxes indicate application of the quadratic Heun- 
to-Heun transformation ( |2.7[ ). If this transformation can be applied after an inde- 
composable Gauss-to-Heun transformation, the relevant parameter specializations 
and composite coverings are indicated to the right (or near the lower right corner) 
of the respective box. If (2.7 1 can be applied after a composite Gauss-to-Heun 
transformation, this is indicated by the => symbol to the right of the H number of 
the composite covering (and to the right of the respective vertical line). 

Some boxes of the same kind touch each other, but that does not have a particu- 
lar meaning. The box for the quadratic E — > HE transformation (2.8) is connected 
to all E — ► E boxes, since this transformation can always be applied without restric- 
tions on the exponent differences. The box for the quadratic E — > E transformation 
(2.6) is connected to all E — > HE boxes, except for the isolated three. 

To show completeness of Figure [3J one must: 
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o Consider all transformations of Tables [2j [3] to which the quadratic Heun-to- 
Heun transformation (2.7) can be applied; and after computing and examin- 
ing the resulting compositions, keep only indecomposable transformations. 

o Check the classical E — > E transformations of [5^1 Table 1] and the qua- 
dratic E — > E transformations, to determine whether the pulled-back hyper- 
geometric equation can ever have exponent differences of the form 1/fc, 1/1, 
consistent with (3.7); if so, composition with a one-parameter E — > HE 
transformation of Table [2] or [3] may be possible. 

o If a pair of E — > E and E — > HE boxes is not connected by a vertical line, 
check that the respective transformations cannot be composed. 

o Check completeness of coverings for each vertical line. 

o Check possible compositions with the Heun-to-Heun transformations of de- 
grees 2 and 4. 

The information of Figure [3] is given in the rightmost columns of Tables [T] [2] [3] 
The compositions are spelled out more explicitly in [301 Appendix B]. A multiple 
occurrence of a covering in Figure [3] means either that it can be decomposed in 
more than one way (as for #3, #5, Hq, #19, #31, #39, #41); or that it appears in 
more than one composition (as for #25, H 2 s, #35, H 3S ); or both (as for #20)- 

The following cases are worth attention. Firstly, there are three ways to de- 
compose the quartic covering #31 in Table [TJ on account of its special branching 
pattern 2+2 = 2+2 = 2+2: 

E(a, a, 2/3) A HE(2a, 2a, 2/3, 2/3), 

(4.3) #31 : <( £(1/2, a, /3) E(2a, /3, /3) A HE{2a, 2a, 2/3, 2/3), 

#£(1/2,1/2, 2a, 2/3) ^- #£(2a, 2a, 2/3, 2/3). 

This is indicated by the 2 x 2 in the rightmost column. The covering #31 occurs 
as a part of the larger compositions #5 and #41 ; see their composition lattices in 
[3U1 (B.5), (B.4)]. Besides, the covering #31 induces the degree 4 Heun-to-Heun 
transformation (4.2). 

The transformation £7(1/2,1/4, a) <A- #£(1/2,1/2, 2a, 2a) is induced by two 
distinct coverings: #31 and #35. Induced by #31, this transformation is the /3 = 1/4 
specialization of (4.3); induced by #35, this transformation is a new one suggested 
by the branching pattern given in Table [3] Both transformations have the factor- 
ization 

(4.4) £(1/2, 1/4, a) A #£(1/2, 1/2, 1/2, 2a) #£(1/2, 1/2, 2a, 2a), 

but they have different sets of t parameters. Both # 31 and #35 appear as parts of 
the degree 8 composite transformation #41. 




5. Existence and uniqueness of coverings 

This section presents an elegant way to conclude that there are no Belyi coverings 
with some branching patterns. The idea is to deduce a pull-back transformation of 
Fuchsian equations that is not possible, because it would relate an equation with 
finite monodromy to an equation with infinite monodromy group, or the pulled- 
back equation would not exist. We apply this idea to all cases of non-existent 
coverings of Tables [T] [2j [3j Moreover, in j ]5.3| this approach is applied to most cases 
of non-existent coverings in the Miranda-Persson list [T5] of K3 elliptic surfaces. 
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As an immediate example, consider the non-existent covering of degree 4 in 
Table [I] If it would exist, the specialization a = 1/2 would give a pull-back from 
E(l/2, 1/2, /?) to a Fuchsian equation with two singularities and (generally) non- 



equal exponent differences ,3,3/3 at them, contradicting part (ii) of Lemma 5.1 
below. Or one can further specialize f3 = 1 or (3 = 1/3 and get a contradiction 
with part (i) of the same lemma. In j]5.1|we prove several assertions from which we 



make non-existence conclusions. Table |5] outlines the non-existence proofs. In § |5.4| 
we seek to show uniqueness (up to Mobius transformations) of the Belyi coverings 
with the encountered branching patterns, by considering implied pull-backs between 
Fuchsian equations with finite monodromy groups. 

5.1. Principal lemmas. The easiest way to conclude non-existence of a Belyi 
covering with a certain branching pattern is to deduce a pull-back transformation 
to a non-existent Fuchsian equation. Here are two basic situations. 

Lemma 5.1. (a) There is no Fuchsian equation on P 1 that has exactly one 
relevant singular point. 
(6) If a Fuchsian equation on P 1 has exactly 2 singular points, their exponent 
differences are equal. 

Proof. If a Fuchsian equation has just one relevant singularity, we can move it to 
infinity and make all points in C ordinary. The differential equation then has the 
form y" + Py' + Qy = 0, where P, Q are polynomials (in the differentiation variable 
x). If P = Q = 0, then the local exponents at the infinity are 0, —1, thus x = oo 
will be an irrelevant singularity. Otherwise x = oo is an irregular singularity, and 
the equation will not be Fuchsian. 

If a Fuchsian equation has 2 singularities, we can assume them to be x = 0, 
x = oo. The Liouville normal form of the equation is then x 2 y" = cy with c G C. 
The exponent differences of this equation equal yT + 4c at both singular points. □ 

Another type of non-existent transformation is a pull-back of a hypergeometric 
equation with finite monodromy to a hypergeometric equation with infinite mon- 
odromy. (A Fuchsian equation has finite monodromy if and only if its solution space 
has a basis consisting of algebraic functions.) The following lemma characterizes 
some hypergeometric equations with finite (or infinite) monodromy groups. 

Lemma 5.2. Consider a hypergeometric equation E = E(a,f3,"f) on P 1 . 

(a) Suppose that a,/3,7 are rational numbers, each having denominator 3. 
Then the monodromy of E will be finite if and only if the sum of the nu- 
merators of a, /3, 7 is even. 

(b) If a is a half-odd-integer, and /3,7 are rational numbers, each having de- 
nominator 4, then the monodromy of E is not finite. 

(c) Suppose that a,/3,7 are integers. Then the monodromy of E will be trivial 
if and only if the sum a+/3+"f is odd, and the triangle inequalities 7 < tt+f3, 
(3 < a + 7, a</3 + 7 are satisfied; otherwise the monodromy is not finite. 

(d) Suppose that a is an integer while /3,7 are half- odd-integers. The set 
{\(3 — 7I ,/3+7} contains two integers of different parity; let k be the integer 
in this set such that k + a is odd. Then the monodromy group of E will 
be isomorphic to Z/2Z if and only if k < a; otherwise the monodromy will 
not be finite. 
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Proof. We use the Schwarz classification of hypergeometric equations with finite 
monodromy for the first two statements; see [33] or §2.7.2]. The only possible 
projective monodromy in statement (a) is the tetrahedral group A4. There are two 
Schwarz types (II and III) for this group: the hypergeometric equation must be 
contiguous either to £7(1/2, 1/3, 1/3) or to £7(1/3, 1/3, 2/3). We must have the latter 
Schwarz type III. Contiguous relations shift the exponent differences by integers 
whose sum is even. That does not change the parity of the numerator sum (of the 
three integers divided by 3), even if an exponent difference is multiplied by —1. 

We do not find the denominator pattern of the statement (6) in the Schwarz list. 
In particular, the two Schwarz types (IV and V) for the octahedral group £4 are 
contiguous to £(1/2, 1/3, 1/4) or £(2/3, 1/4, 1/4). 

For the claim (c), a representative solution of the generic hypergeometric equa- 
tion with trivial monodromy is 2Fi( — n, £+1; —n — m \ z), with n,m,i non-negative 
integers; see [27J Theorem 2.4(5)]. Up to a permutation, one has that a = n+m+1, 
/3 = n + £+l, j = m + e + l; that is 

_ q + g-7-1 _ a+j-p-1 (3 + j-a-l 

{o.i) n — - , m — ^ ' 2 

If one of these three numbers is a negative integer, the singular point with the largest 
exponent difference is logarithmic [37J § 9] . If each of the above three numbers is a 
half-odd- integer, all three singular points are logarithmic [37J §5]. 

The assertion (d) is a reformulation of [2"51 Theorem 5.1], stated in the context of 
hypergeometric equations with either logarithmic solutions or the Z /2Z monodromy 
group. □ 

Existence (and uniqueness) of coverings with a given branching pattern can also 
be decided on the basis of transformations of some hypergeometric equations with 
infinite monodromies. The following lemma implies that there are no transforma- 
tions of £(1/2,1/4,1/4) into itself of degrees 6, 12, 14, 21, 22, 24, or generally, 
of degrees = 3 (mod 4), even if suitable branching patterns of these degrees exist. 
Similarly, there are no transformations of £7(1/2, 1/3, 1/6) or £(1/3, 1/3, 1/3) into 
themselves of degrees 6, 10, 15, 18, 22, 24, or generally, of degrees = 2 (mod 3). 
This lemma eludicates the non-uniqueness of H44 and i?2i- 

Lemma 5.3. (a) Up to Mobius transformations, the number of degree- D pull- 
back coverings of £7(1/2, 1/4, 1/4) into itself is equal to the number of integer 
solutions (a, b) with a ^ 0, b > 0, of the equation D = a 2 + b 2 . 
(b) Up to Mobius transformations, the number of degree- D pull-back coverings 
of £7(1/2, 1/3, 1/6) or £7(1/3, 1/3, 1/3) into itself is equal to the number of 
integer solutions (a, b) with a ^ 0, b > a, of the equation D = a 2 — ab + b 2 . 

Proof. According to [3§1 §8], the transformations of £7(1/2,1/4,1/4) into itself 
correspond to isogenics of the j = 1728 elliptic curve y 2 — x 3 — x. The ring of 
isogenies is isomorphic to the ring Z[i] of Gaussian integers, and the degree of a pull- 
back is equal to the norm a 2 +b 2 of the corresponding a+bi. In particular, the trivial 
and fractional- linear transformations correspond to the units ±1, ±i. Therefore one 
must count a + bi & Z[i] such that \a + bi\ 2 = D and arg(a + hi) G [0, tt/2). 

Similarly §8], the transformations of £(1/2, 1/3, 1/6) or £7(1/3, 1/3, 1/3) into 
itself correspond to isogenies of the j = elliptic curves y 2 — x 3 — 1 or x 3 + y 3 = 1. 
The ring of isogenies is isomorphic to the ring of Eisenstein integers Z [u] . The degree 
of a pull-back is equal to the norm a 2 ~ab + b 2 of the corresponding a + but. Trivial 
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Table 5. Unrealizable branching patterns, with a proof indication 



or Mobius transformations correspond to the units ±l,±u;,±(aj + l). Therefore one 
must count a + buj e Z[w] such that \a + bu\ — D and arg(a + bco) S [0, tt/3). □ 

5.2. Nonexistence of coverings. Tables [2] [3] have 27 entries with nonexistent 
Belyi coverings. One branching pattern appears twice among the type (2,4) candi- 
dates, hence the two tables actually have 26 different branching patterns with no 
covering. They are labeled Ni, . . . , N26. The repeating branching pattern is labelled 



/V23. Nonexistence is in each case an immediate consequence of some lemma in 1 5.1 
Mostly by specialization of the free parameter, one either derives a pull-back from 
a hypergeometric equation to a nonexistent Fuchsian equation, or a pull-back of a 
hypergeometric equation with finite monodromy to a hypergeometric equation with 
infinite monodromy, or a nonexistent pull-back of E(l/2, 1/3, 1/6) into itself. The 
unrealizable branching patterns and the applicable lemmas are listed in Table [5j 

The non-existent covering of Table [T] is given the last number /V27. Its non- 
existence was already demonstrated at the beginning of this section. 

Only for N 2 i and N 2 s the used implied transformation is not a specialization of 
a respective Gauss-to-Heun pull-back of the classification in pi To prove N 2 i by 
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the specialization a = 1/5, one would need to inspect the 10 icosahedral Schwarz 
types in § 2.7.2]. The case N 2 3 can be proved using the specialization a = 1/4 of 
either of the two candidate transformations in Table [3j by invoking Lemma |5.3[ b). 
Note that to use a hypergeometric equation with only two relevant singularities, 
one must ensure that it is of the form E(l, a, a). In particular, Lemma |5.1| (&) 
does not apply to the branching covering [2]6=[3]4=9+l+l+l and its pull-backs 
from £'(1/2,1/3,1), because logarithmic singularities rather than ordinary points 
appear. And indeed, the covering H± exists. 



5.3. The Miranda-Persson classification. The lemmas of j |5.1| can be applied 
to the problem of the existence of Belyi maps that would yield semi-stable elliptic 
fibrations of K3 surfaces with 6 singular fibers, sorted out by Miranda, Person [18] 
and Beukers, Montanus ; 4j . The degree of the relevant Belyi maps is 24, and their 
branching patterns have the form [2] 12 = [3]§ = P, where P = a+b+c+d+e+f is 
a partition of 24 with exactly 6 parts. There are 199 of these branching patterns 
in total. Miranda and Persson [TS] proved that Belyi coverings (and elliptic K3 
surfaces) exist in 112 cases, and do not exist in the remaining 87 cases. Beukers 
and Montanus [4] computed alQ those Belyi maps and checked non-existence for 
the 87 partitions. 

The non-existence proof in [TS] broadly relies on two techniques. First, Miranda 
and Persson widen the space of considered branching patterns to include partitions 
P with more than six part^Jand conclude non-existence of coverings for a partition 
a\ + . . . + a s from non-existence for a partition <zi + . . . + a s —i + a' s + a" with 
a s = a' s + a", using 18, Lemma (2.4)]. Secondly, they get contradicting conclusions 
about the torsion of the assumed elliptic surfaces in several non-existing cases. 
In [3], non-existence is concluded either by using a sum over the characters of ^24 
that counts coverings (not necessarily connected, with some rational weights) with 
a given branching pattern, or by direct computation. Let E denote the counting 
character sum just mentioned, given in [4j Theorem 3.2]. The large table in [4] does 
not list the 47 partitions (out of the total 87) for which E = 0. 

Here we show that most of the non-existent cases in the Miranda-Persson list 



can be deduced using the methods of [5.1 Here are 22 partitions out of the 40 



ones with E ^ for which the non-existence can be proved by using Lemmas 5.1 
[5~2] |Q] directly: 

14+[2] 5 , 9+[3] 5 , 15+[2] 4 +l, 13+3+[2] 4 , 12+4+[2] 4 , ll+5+[2] 4 , 10+6+[2] 4 , 



^As pointed out in the AMS MathSciNet review by David P. Roberts, the table in [4] omits 
one Belyi covering for the partition 10+6+4+2+1-1-1. Our computation confirms existence of two 
(rather than one) Belyi coverings for this partition: 

(144a: 8 + 384a: 7 + 1120a: 6 - 784a: 3 + 756a: 2 - 240a; + 25) 3 

108a: 6 (14x - 5) 4 (Ax - l) 2 (9a; 2 + 24a; + 70) ' 
(144a; 8 - 1536a; 7 + 5248a; 6 - 5568a; 5 - 720a; 4 + 512a; 3 + 192a; 2 + 24a; + l) 3 
108 (82 + l) 6 a; 4 (x - 3) 2 (9a; 2 - 42a; - 5) 

The second covering is missing in the Beukers— Montanus list. In total, there are 59 branching 
patterns (among the 112 indicated by Miranda and Persson) with a unique Belyi map up to 
Mobius transformations; 125 Galois orbits of the Belyi maps, of size at most 4; and 191 different 
Belyi maps or dessins denfant. 

2 Therefore coverings with more than 3 branching fibers are allowed. Instead of the coverings, 
permutation representations of their monodromy are considered in I18| . 
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ll+[3] 4 +l, 10+[3] 4 +2, 8+4+[3] 4 , 13+4+[2] 8 +l, ll+6+[2] 3 +l, ll+4+3+[2] 3 , 
10+4+4+[2] 3 , 9+8+[2] 3 +l, 9+6+3+[2] 3 , 8+7+3+[2] 3 , 8+5+5+[2] 3 , 7+7+4+[2] 3 , 
10+[4] a +l+l, 6+[4] 3 +3+3, [6] 3 +3+2+l. 

The choice of the starting £7(1/2, 1/3, 1/k) that yields a non-existent covering is 
indicated by the [k) n notation. Next, here are 22 partitions out of the 47 ones with 
£ = to which our lemmas apply directly: 

9+7+[2] 4 , 7+5+[3] 4 , 7+[4] 4 +l, 6+[4] 4 +2, 5+[4] 4 +3, [5] 4 +3+l, 
9+5+4+[2] 3 , 7+6+5+[2] 3 , 13+[3] 3 +l+l, ll+[3] 3 +2+2, 10+4+[3] 3 +l, 
8+6+[3] 3 +l, 8+5+[3] 3 +2, 7+7+[3] 3 +l, 7+6+[3] 3 +2, 7+4+4+[3] 3 , 6+5+4+[3] 3 , 
5+5+5+[3] 3 , 9+[4] 3 +2+l, 6+5+[4] 3 +l, 7+[4] 3 +3+2, 5+5+[4] 3 +2. 

Additionally, the four cases 7+[5] 3 +l+l, 6+[5] 3 +2+l, [5] 3 +4+4+l, [5] 3 +4+3+2 
with £ = are concluded by inspecting the icosahedral hypergeometric equations 
in the Schwarz table O § 2.7.2]. In total, this shows 48 out of the 87 cases. 

More cases of non-existence can be deduced from implied pull-backs to Fuchsian 
equations with 3 non-apparent singularities and a few apparent singularities. These 
equations are gauge "contiguous" to target hypergeometric equations (with infinite 
or infinite monodromy) as the local exponent differences differ at all points by inte- 
gers. The total shift of the exponent differences, including those from the difference 
1 for ordinary points of hypergeometric equations, must be an even integer. In this 
way, non-existence for the following 7 partitions with E / can be shown: 

10+6+[3] 2 +l+l, 9+9+[3]i + l+l+l, 8+6+[3] 2 +2+2, 7+6+6+[3]i+l+l, 
7+6+4+[3] 2 +l, 6+5+5+[3] 2 +2, 8+6+[4] 2 +l+l. 

In each case, the apparent singularities are represented by the branching orders 
that are integer multiples of the bracketed numbers. And here are 7 partitions with 
£ = that can be handled in the same way: 

9+7+[3]a+l+l, 9+5+[3] 2 +2+2, 9+4+4+[3] 2 +l, 6+6+5+[3]i+2+2, 
6+6+4+4+[3]i+l, 8+5+[4] a +2+l, 8+[4] 2 +3+3+2. 

Besides, a pull-back from £(1/2, 1/3, 1/3) can be applied to show the non-existence 
for 9+6+6+1+1+1, with E It is trickier to combine parts (c), (d) of Lemma 



5.2 with gauge shifts. 



Of the remaining 87 — 48 — 7 — 7—1 = 24 partitions, the following 6 (with E / 0) 
and 11 (with £ = 0) partitions could be handled with a full knowledge of Heun 
equations with finite monodromy (that are not classified yet): 

10+8+[2] 2 +l+l, 13+[4] 2 +l+l+l, ll+[4] 2 +2+2+l, 9+[4] 2 +3+2+2, 
9+[5] 2 +3+l+l, 8+[5] 2 +4+l+l; 
9+6+4+[2] 2 +l, 8+8+3+[2] 2 +l, 7+7+5+[2] 2 +l, 7+6+4+3+[2] 2 , 
10+5+[3] 2 +2+l, 7+5+5+[3] 2 +l, 9+5+[4] 2 +l+l, 7+5+[4] 2 +2+2, 
8+[5] 2 +3+2+l, 6+[5] 2 +4+3+l, 6+[5] 2 +4+2+2. 

Besides, a pull-back from E(l/2, 1/3, 1/4) could be then applied to two partitions 
with £ = 0: 12+8+1+1+1+1, 8+8+5+1+1+1. Other 3 partitions (with E/0) 

12+5+[4]i+l+l+l, 10+[5]i+4+3+l+l, 9+8+[4]i+l+l+l, 

could be decided by Fuchsian equations with 4+1 singularities (i.e., 4 non-apparent 
and 1 apparent). There remain only two partitions: 7+7+6+2+1+1 with E / 0, 
and 7+7+4+3+2+1 with £ = 0. Their non-existence might be decided by using im- 
plied pull-backs from E(l/2, 1/3, 1/2) to Fuchsian equations with 4+1 singularities 
and the monodromy group D 2 or Z/2Z. 
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5.4. Uniqueness of coverings. Uniqueness of Gauss-to-Heun transformations 
(and of their coverings) with a plausible branching pattern can be concluded from 
uniqueness of specialized Gauss-to-Gauss transformations. In particular, the cov- 
erings H^-H^b, H47 appear in the classical hypergeometric transformations 
listed by Goursat [8]. The coverings Hi, H2, H?, H$, Hn, His, H42 appear in the 
hypergeometric transformations from E(k, £, m) with k, £, m positive integers satis- 
fying l/k + l/£+l/m < 1. As determined in (and §9]), these pull-backs are 
unique up to Mobius transformations as well. The coverings H31, H39, H41, 
apply to hypergeometric transformations from E(l/2, 1/2, a) with infinite dihedral 
monodromy [311 §4]. The pulled-back equations have infinite cyclic or dihedral 
monodromy. They are, respectively, 

£(l,2a,2a), #(1/2, 1/2, 6a), £(l,4a,4a), E(l/2, 1/2, 5a). 

The cyclic covering gives the pull-back E(l,a,a) J?(l,4a,4a) of hyper- 
geometric equations with infinite cyclic monodromy. 

Non-unique Gauss-to-Gauss transformations appear when hypergeometric equa- 
tions E(k,£,m) are pulled-back, with k,£,m positive integers satisfying l/k + l/£+ 
1/m ^ 1. It the equality holds, these hypergeometric functions are integrals of 
holomorphic differentials on j = 1728 or j = elliptic curves [29j §8]. Lemma 



5.3 counts the coverings H 3 , i? 12 , H 2 \, H i0 , if 44 , H ie . If can be established (by 
identifying transformations of holomorphic differentials on the curves y 2 = x 3 — 1 
and x 3 + y 3 = 1, y 2 = x 6 + 1) that the transformations from E(l/2, 1/3, 1/6) to 
-E(l/3, 1/3, 1/3) or E(2/3, 1/6, 1/6) are compositions of the pull-backs of Lemma 
|5.3| with quadratic transformations. This applies to the coverings -H15, -H19, i?38- 

The hypergeometric equations E(k, £, m) with 1/k + l/£ + 1/m > 1 have finite 
monodromy groups. The hypergeometric solutions are thereby algebraic functions. 
These equations play a fundamental role in the classical heory of algebraic solutions 
of second order Fuchsian equations: 

o E(l, 1/k, 1/k), with the finite cyclic monodromy Ck- 

o £^(1/2, 1/2, 1/k), with the dihedral projective monodromy Dk- 

o _E(l/2, 1/3, 1/3), with the tetrahedral projective monodromy At. 

o E(l/2, 1/3, 1/4), with the octahedral projective monodromy S4. 

o E(l/2, 1/3, 1/5), with the icosahedral projective monodromy A$. 

By a celebrated theorem of Klein [TJ, all second order Fuchsian equations on P 1 
with a finite monodromy group are pull-backs of one of these standard hypergeo- 
metric equations, with the same projective monodromy group. These Klein trans- 
formations are known to be unique up to Mobius transformations [1 . However, pull- 
back transformations between hypergeometric equations with different projective 
monodromy need not to be unique. Litcanu |151 Theorem 2.1] noted non-uniqueness 
of the pull-backs from £(1/2,1/3,1/4) to £(1/2,1/2,1/2) and £(1,1/2,1/2), of 
degree 6 and 12 respectively. The non-uniqueness is caused by pairs of different 
branching patterns though, e.g., 2+2+2=3+3=2+2+2 and 2+2+1+1=3+3=4+2. 
The example of £(1/2,1/2,1/5) £(1/2,1/2,2) in [3TJ §5.4] shows that non- 
unique coverings with the same branching pattern easily occur for pull-backs to 
equations with apparent singularities. Besides, many compositions of 

(5.2) H 37 : £(1/2, 1/3, 1/5) A £(1/2, 1/3, 1/3) 
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with transformations from the tetrahedral equation are not unique either, because 
the properly normalized £37 is defined over Q(\/— 15); see formula [221 (50)]. 

In Table[4j the coverings £9, £io, £13, £16, £22, £24 give Klein transformations 
of £7(1/2, 1/3, 1/5) to the following hypergeometric equations, respectively: 

£7(1/3,1/5,4/5), £7(1/3,2/5,3/5), £7(1/2,1/5,3/5), 
£(2/3,1/5,2/5), £(1/2,1/3,2/5), £(1/3,2/3,1/5). 

This illustrates the Schwarz types VIII, XV, IX, X, XIV, XII, respectively. The 
other icosahedral Schwarz types are represented by £(1/3, 1/3, 2/5), £(1/5, 1/5, 4/5), 
£(2/5,2/5,2/5), and the standard £(1/2,1/3,1/5). Uniqueness of the coverings 
£14, £17, £23, £29, £30 is established by noting these Klein transformations: 

£(1/2, 1/3, 1/3) A £(1/2, 2/3, 4/3), £(1/2, 1/2, 1/3) A £(3/2, 3/2, 2/3), 

£(1/2, 1/3, 1/4) A £(1/2, 1/3, 3/4), £(1/2, 1/3, 1/4) A £(1/2, 2/3, 1/4), 

£(1/2, 1/3, 1/3) A £(1/2, 2/3, 2/3). 

These considerations of reduction to hypergeometric transformations do not im- 
mediately establish uniqueness of 10 coverings in Table [4] Those coverings induce 
rather attractive transformations between hypergeometric equations with different 
finite monodromy. In particular, £g, £28 pull-back £(1/2,1/3,1/3) to £(1,1,1) 
and £(1,1/2,1/2); then £ 5 , £ 20 , £25, £27, £36 transform £(1/2,1/3,1/4) to 
£(1,1/2,1/2), £(1,1/3,1/3), £(1/2,1/3,2/3), £(1/2,1/2,1/2), £(1/2,1/2,1/3), 
respectively; and finally, £4, £ 26 , £37 pull-back £(1/2, 1/3, 1/5) to £(1, 1/5, 1/5), 
£(1/2,1/2,1/5) and |5T2| . Many of the coverings pull-back £(1/2,1/3,1/2) or 
other dihedral hypergeometric equations to hypergeometric equations with simpler 
dihedral or cyclic monodromy. 
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